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Abstract 

This  paper  presents  investigations  on  the  dynamic  behaviour  of  DMFCs.  Experimental  cell  voltage  responses  to  cell  current  step  changes 
are  investigated  with  transfer  function  analysis.  The  voltage  response  in  most  places  shows  significant  overshooting  behaviour  (up  to  80%) 
when  applying  cell  current  steps  to  the  DMFC.  Linear  system  analysis  is  applied  to  a  set  of  dynamic  models  with  different  complexity,  and  a 
criterion  for  cell  voltage  overshooting  is  developed.  In  all  models  the  dominating  physico-chemical  phenomenon  is  the  transport  of  methanol 
through  the  membrane,  causing  cathodic  overpotential  overshooting.  This  is  supported  by  experimental  analysis.  The  proposed  models  and 
transfer  functions  are  suitable  for  controller  design. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

First  direct  methanol  fuel  cell  (DMFC)  systems  are  com¬ 
mercially  available  (e.g.  from  SmartFuelCell,  Germany)  after 
more  than  one  decade  of  research.  They  are  mainly  designed 
for  the  use  as  mobile  supplies  of  electrical  power  in  the  range 
of  up  to  a  few  hundred  Watts.  Nonetheless,  many  problems 
related  to  the  operation  of  DMFCs  under  realistic,  i.e.  dy¬ 
namic,  operating  conditions  have  yet  to  be  solved.  Therefore, 
all  available  fuel  cell  systems  currently  use  a  large  battery  or 
supercapacitor  system  between  electric  load  and  DMFC  to 
buffer  dynamic  load  demands.  Usually,  the  DMFC  is  oper¬ 
ated  at  a  moderate,  constant  electric  load  which  is  suitable 
for  long-term  performance.  Thus,  the  DMFC  provides  for  a 
high  state  of  charge  of  the  much  larger  battery  system.  To 
widen  the  market  for  the  DMFC,  though,  further  miniaturi¬ 
sation  and  lower  system  costs  have  to  be  achieved.  Therefore, 
it  is  attractive  to  be  able  to  operate  the  DMFC  system  reliably 
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and  effectively  without  such  a  buffer  system,  or  at  least  with  a 
much  smaller  one.  This  in  turn  requires  a  deeper  understand¬ 
ing  of  the  behaviour  of  the  DMFC  under  realistic  dynamic 
operating  conditions,  like  changes  in  the  cell  current  (i.e.  the 
external  load  demand)  and  the  methanol  feed  concentration 
(as  one  of  the  few  controllable  system  variables).  Reduced 
but  realistic  mathematical  models  can  prove  useful  in  un¬ 
derstanding  dynamic  interactions  between  the  main  physical 
and  chemical  phenomena,  and  hence  may  act  as  a  basis  for 
improved  cell,  system  and  controller  design. 

In  recent  years,  research  has  been  conducted  on  experi¬ 
mental  DMFC  responses  to  dynamic  changes  in  cell  current 
[1^4].  In  the  two  latter  publications,  the  authors  performed 
various  cell  current  steps  and  waited  for  a  few  seconds  until 
performing  the  next  step.  They  also  conducted  gradual  cell 
loading  and  unloading  from  or  to  OC  V.  In  addition,  they  com¬ 
bined  cell  current  steps  and  gradual  loading  and  unloading 
to  driving  cycles.  Such  measurements  proved  useful  for  em¬ 
pirical  one-step-ahead  prediction  models  of  DMFC  stacks 
[2].  Nonetheless,  the  experiments  are  not  ideal  for  analy¬ 
sis  of  processes  occurring  inside  the  fuel  cell.  For  example, 
the  observation  time  after  a  current  step  is  mostly  only  30  s. 
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Nomenclature 

A 

system  matrix 

As 

geometric  electrode  area  (26  x  10-4  m2) 

b\ ,  b2 

coefficients  for  the  Laplace  transformed  flux 

B 

input  matrix 

c 

concentration  (mol  m-3) 

cref 

lCH3OH 

reference  methanol  concentration 
(1000  mol  m-3) 

JA 

CH+ 

proton  concentration  in  membrane  pores 
(1200  mol  m-3)  [7] 

C 

output  matrix 

cAC 

anode  double  layer  capacity  (1827 Cm-2) 
(own  CV  measurements) 

/'-•AC 

cCHsOH 

concentration  variable  in  Laplace  domain 

Ccc 

cathode  double  layer  capacity  (907  C  m-2) 
(own  CV  measurements) 

D 

direct  transmission  matrix 

dAC 

thickness  of  anode  catalyst  layer  (35  x 
10-6  m)  [10] 

dAD 

thickness  of  anode  diffusion  layer  (1.7  x 
10-4m)  [10] 

dM 

thickness  of  (fully  hydrated)  Nafion™  105 
membrane  (10-4  m)  [10] 

/A(i  [,  oi  i 

diffusion  coefficient  of  methanol  in  water  at 
333  K  (3.187  x  10-9m2s-1)  (calc.  acc.  [24], 

p.  600) 

^CHsOH 

diffusion  coefficient  of  methanol  in  membrane 
at  333  K  (4.7  x  10-1°m2s-1)  (own 
measurements) 

dm 

u\i+ 

diffusion  coefficient  of  protons  in  membrane 
(4.5  x  10-9  m2  s-1)  [25] 

F 

Faraday  constant  (96485  C  mol-1) 

Fa 

volume  flow  rate  of  anode  feed 
(0.5  dm3  min-1) 

Fi,F2  ( F]  rat'  ^2. rat)  transcendental  parts  of  the 

membrane  transfer  function  (rational  forms) 

G 

transfer  function 

i 

counter 

*cell 

cell  current  density  (A  m-2) 

*max 

maximum  initial  cell  current  density  for 
overshooting  behaviour  (A  m-2) 

/ 

current  variable  in  Laplace  domain 

j 

imaginary  unit 

k 

constant  in  transfer  function 

kAD 

effective  mass  transport  coefficient  in 

anode  diffusion  layer  (  Dch3 oh (<? A D ) 1  "V<7 AD ) 

(ms-1) 

k(p 

electrokinetic  permeability  of  membrane 
(1.13  x  10-19m2)  [25] 

k<l>  2 

electro-osmotic  parameter 
(1.33  x  10-3  m3  C-1)  (see  text) 

m 

number  of  time-relevant  singularities 

n 

number  of  poles  and  zeros 

»ch3oh 

flux  density  of  methanol  (mol  m  2  s  1 ) 

Z'max.Pt 

number  of  COx  adsorption  sites  on  Pt 
(0. 1 1  mol  CO  m-2) 

/Vcihou  variable  of  methanol  flow  in  Laplace  domain 

P 

pole  of  a  transfer  function 

Pe 

Peclet  number  (k02£/Mz'ceii/-DcH3OH) 

^ai(^aio)  reaction  rate  (constant)  for  anodic  methanol 

oxidation  (mol  m  2  s  1 ) 

rA2(r,\2o)  reaction  rate  (constant)  for  anodic  COx 

oxidation  (mol  m-2  s- 1 ) 

n 

n 

o 

reaction  rate  (constant)  for  cathodic  reaction 
(0.5  x  10-4molm-2  s-1) 

R 

universal  gas  constant  (8.314  Jmol-1  K-1) 

Re  1 

Ohmic  resistance  of  the  membrane  (<7m/zfm) 
(£2m2) 

res 

residue 

s 

Laplace  variable 

t 

time  (s) 

T 

cell  temperature  (333  K) 

u(U) 

input  variable  (in  Laplace  domain) 

TjO 

standard  cell  voltage  (1.213  V) 

(thermodynamic  calculation) 

^cell 

cell  voltage  (V) 

yAC 

volume  of  anode  catalyst  layer 
(7.37  x  10-8m3) 

wAC 

anode  (and  cathode)  catalyst  loading 
(5  mg  cm-2) 

x(X) 

state  variable  (in  Laplace  domain) 

y(Y) 

output  variable  (in  Laplace  domain) 

z 

spatial  coordinate  in  membrane  (m) 

Z 

zero  of  a  transfer  function 

Greek  letters 

(j,\  i ,  a A2  charge  transfer  coefficient  for  anodic  reactions 

(0.5) 

ac 

charge  transfer  coefficient  for  cathodic  reaction 
(0.5) 

f-AC 

porosity  of  anode  catalyst  layer  (0.81)  [10] 

fAD 

porosity  of  anode  diffusion  layer  (0.71)  [10] 

porosity  of  membrane  (0.4)  [10] 

V 

electrode  overpotential  (V) 

km 

ionic  conductivity  of  membrane 

nief 

&COx 

(0.0873  f2-1  m-1)  (calc.  acc.  [12]) 
reference  COx  surface  coverage  of  active  sites 
of  Pt  catalyst  (0.5) 

Ocox 

COx  surface  coverage  of  active  sites  of 

Pt  catalyst 

M 

pore  fluid  viscosity  in  membrane 
(4.823 kgm-1  s-1)  (calc.  acc.  [24], 
p.  441,  455) 

P 

radius  in  complex  plane 

Subscripts 

A 

anode 

AC 

anode  catalyst  layer 

U.  Krewer,  K.  Sundmacher  /  Journal  of  Power  Sources  154  (2005)  153-170 


155 


AE 

anodic  overpotential 

AT 

surface  coverage 

C 

cathode 

CH3OH  methanol 

COx 

adsorbed  intermediates  occurring  during 
methanol  oxidation,  e.g.  CO,  COH,  CHO 

num 

numerically  evaluated 

rel 

relative 

ss 

at  steady  state  conditions  prior  to  current  step 

tot 

total 

Superscripts 

A 

anode  compartment 

AC 

anode  catalyst  layer 

AD 

anode  diffusion  layer 

CC 

cathode  catalyst  layer 

M 

membrane 

M,  A/(M,  C)  boundary  membrane-anode/ 

(membrane-cathode) 

This  time-span  is  sufficient  for  the  above  mentioned  one-step- 
ahead  models,  but  it  is  not  sufficient  for  reaching  a  steady  state 
for  all  dynamic  variables. 

This  publication  presents  a  systematic  approach  towards 
recording  and  analysing  the  DMFC  voltage  responses  to  load 
steps  between  two  steady  states.  First,  a  polarisation  curve 
is  recorded,  and  the  characteristic  operating  regimes  (e.g. 
pseudo-ohmic  regime)  are  identified.  This  is  followed  by  var¬ 
ious  current  steps  between  the  identified  regimes.  The  voltage 
response  is  recorded  for  600  s.  A  subsequent  analysis  shows 
the  characteristic  influence  of  the  operating  regimes  on  the 
cell  response. 

Furthermore,  mathematical  models  of  different  complex¬ 
ity  are  formulated  and  evaluated  using  the  experimental  data. 
A  variety  of  models  for  different  fuel  cell  types  exists  in  lit¬ 
erature,  most  are  steady  state  models.  Dynamic  models  are 
published  e.g.  by  Heidebrecht  and  Sundmacher  [5]  (MCFC), 
Ceraolo  et  al.  [6]  (PEMFC),  Sundmacher  et  al.  [7]  (DMFC), 
and  Yerramalla  et  al.  [8]  (PEMFC). 

Only  few  groups  investigate  on  the  dynamics  of  the 
DMFC.  Besides  empirical  dynamic  DMFC  models,  e.g. 
[2],  the  dynamic  behaviour  of  DMFCs  to  concentration 
changes  was  modelled  [7,9].  Zhou  et  al.  [9]  were  among 
the  first  to  apply  system  theoretical  approaches  on  fuel  cell 
models,  and  in  particular  on  a  DMFC  model.  They  showed 
that  non-minimal-phase  behaviour,  which  was  observed  in  a 
concentration  step  experiment,  can  be  predicted  via  analysis 
of  the  transfer  function  of  a  respective  linearised  model. 
Linear  system  analysis  has  many  advantages:  analytical 
analysis  of  models  can  be  conducted  which  may  yield 
principle  understanding  of  the  influence  of  the  underlying 
physico-chemical  phenomena.  Furthermore,  transfer  func¬ 
tions  are  a  basis  for  controller  design.  Controllers  based  on 
physico-chemical  process  models  are  of  particular  interest. 


In  the  following,  dynamic  DMFC  models  will  be  presented 
which  analyse  DMFC  responses  to  dynamic  load  changes 
by  using  classical  methods  of  system  analysis.  It  will  be 
shown  that  the  applied  methods  predict  conditions  for  experi¬ 
mentally  observed  overshooting  behaviour.  Additionally,  the 
influence  of  reaction,  methanol  crossover  and  double  layer 
charging  will  be  investigated. 

The  presented  models  hold  for  analysis  of  an  operating 
regime  in  which  mass  transport  influences  are  minimised: 
low  anodic  residence  time  (i.e.  negligible  methanol  concen¬ 
tration  distribution),  low  current  densities  (i.e.  negligible  CO2 
bubble  influence).  This  allows  a  much  more  indepth  analy¬ 
sis  of  the  processes  occurring  inside  the  MEA.  Future  works 
will  systematically  add  outer  influences  like  mass  transport 
phenomena  and  a  refined  methanol  oxidation  mechanism. 

2.  Experimental  setup 

2.1.  Applied  DMFC  design  and  materials 

The  experiments  were  carried  out  using  a  single  DMFC 
fed  with  air  and  liquid  methanol/water  solutions  [10].  The 
identical  anode  and  cathode  monopolar  plates  are  made  from 
graphite  material  (thickness  7  mm,  material  code  FU4369) 
supplied  by  Schunk  Kohlenstofftechnik  (Germany).  The  nec¬ 
essary  flow  field  structures  for  the  reactant  distribution  over 
the  MEA  surface  are  millcut  into  the  plates.  They  consist  of 
parallel  channels  of  2  mm  width  and  2  mm  depth,  with  1  mm 
wide  ribs  between  them.  A  distributor  and  collector  channel 
connect  the  parallel  channels  to  the  inlet  and  outlet  ports,  re¬ 
spectively.  The  media  (air  and  methanol/water  solution)  are 
supplied  in  one  corner  of  the  rectangular  flow  field  and  leave 
at  the  opposite  corner.  The  flow  field  itself  has  the  outer  di¬ 
mensions  65  mm  x  40  mm,  identical  to  the  catalyst  layer  on 
the  MEAs,  which  leads  to  an  active  area  of  As  =  26  cm2. 

As  diffusion  layer  PTFE-coated  TORAY  carbon  paper 
(TGP-H-060)  is  used,  with  a  PTFE  loading  between  20  and 
25  mass%  with  respect  to  the  uncoated  material. 

Finally,  the  membrane  electrode  assemblies  (MEA)  are 
prepared  from  NAFION™  N-105  membrane  foil,  onto  which 
the  catalyst  layers  are  applied  using  an  airbrush  technique  de¬ 
veloped  by  ZSW,  Ulm  (Germany)  [11].  The  anode  catalyst 
layer  features  a  catalyst  loading  of  5  mg  cm-2  (unsupported) 
platinum  ruthenium  black  (Alfa  Aesar  Johnson  Matthey 
HiSPEC™  6000)  and  a  NAFION™  content  of  15mass% 
relative  to  the  metal  loading  (i.e.  0.75  mg  cm-2).  The  cath¬ 
ode  catalyst  layer  has  the  same  metal  loading,  but  as  catalyst 
(unsupported)  pure  platinum  black  is  used  (Alfa  Aesar  John¬ 
son  Matthey  HiSPEC™  1000)  and  the  NAFION™  content 
is  10  mass%  relative  to  the  metal  loading  (i.e.  0.5  mg  cm-2). 

The  DMFC  is  completed  by  gold-plated  copper  plates 
as  current  collectors  and  stainless  steel  plates  for  bracing 
the  whole  sandwich  structure.  A  torque  of  5  Nm  is  exerted 
on  the  screws,  which  hold  together  the  steel  back  plates. 
After  assembly,  each  DMFC  is  conditioned  and  evaluated 
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by  operation  with  pure  humidified  hydrogen  and  air  for  3  x 
8  h,  before  it  is  operated  on  methanol  solutions. 

2.1.1.  DMFC  miniplant 

For  full-scale  testing  of  DMFCs,  a  miniplant  was  designed 
by  Schultz  [10],  which  is  fully  automated  (process  control 
system  PC-S7AVinCC  by  Siemens)  to  enable  automatic  test¬ 
ing  procedures,  with  a  special  focus  on  dynamic  operation. 
Fig.  1  shows  a  simplified  flowsheet  of  the  miniplant. 

The  DMFC  cathode  is  supplied  with  dry  air  (dew  point 
^271  K)  at  flow  rates  between  0.4  and  S.Oscbmh^1  (mass 
flow  controller  F101,  type  Mass6020  by  Buerkert  AG,  Ger¬ 
many)  at  cathode  outlet  pressures  of  ambient  up  to  5  bars 
absolute  ( 1  (P-5  x  105Pa).  The  air  is  pre-heated  in  a  plate 
heat  exchanger  (W101),  air  temperatures  and  pressures  are 
measured  at  the  cathode  inlet  and  outlet.  At  the  cathode  outlet, 
also  the  relative  humidity  of  the  air  is  measured  (Q202,  type 
HygroClip  IE  by  rotronic  AG,  Switzerland).  Finally  the  cath¬ 
ode  exhaust  air  enters  a  condensor,  where  it  is  dried  to  reach  a 
dew  point  below  263  K  (condensate  is  collected).  The  dry  air 


is  sent  into  a  fume  hood,  while  its  oxygen  and  carbon  dioxide 
contents  are  measured.  The  oxygen  sensor  (Q204)  is  a  para¬ 
magnetic  sensor  (PAROX  1000  H  by  MBE  AG,  Switzerland), 
while  carbon  dioxide  is  measured  using  an  IR-sensor  (Q203, 
type  OEM-NDIR  EGC-5%  by  Pewatron  AG,  Switzerland). 

On  the  DMFC  anode  side,  a  liquid  recycle  loop  is  installed. 
It  consists  of  two  alternative  cycles,  one  for  methanol/water 
solution  and  one  for  pure  water.  The  purpose  of  this  is  to  en¬ 
able  a  pulsed  periodic  operation  of  the  DMFC,  where  the  an¬ 
ode  feed  is  changed  stepwise  between  methanol/water  solu¬ 
tion  and  pure  water  automatically.  Both  branches  of  the  anode 
cycle  feature  vessels  for  pressure  equilibration  and  carbon 
dioxide  removal  (B1  and  B2),  gear  pumps  (P401  and  P402) 
and  heat  exchangers  (W403  and  W406).  Flow  rates  between 
0.3  and  5  dm3  min-1  can  be  achieved.  Automatic  valves 
(V403/V404  and  V408/V409)  enable  a  flexible  and  prac¬ 
tically  immediate  change  between  methanol/water  solution 
and  pure  water  anode  feed  without  causing  significant  distur¬ 
bances  in  liquid  flow  rate  and  pressure.  The  flow  rate  is  mea¬ 
sured  by  a  Coriolis-type  sensor  (F401,  type  MASS  2100  DI6 


i - 1 


Fig.  1.  Simplified  flow  scheme  of  DMFC  miniplant  with  all  important  sensors  and  actuators. 
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by  Danfoss,  Denmark).  Like  on  the  cathode  side,  the  medium 
temperature  and  pressure  are  measured  at  the  cell  inlet  and 
outlet.  The  methanol  concentration  of  the  anode  inlet  medium 
is  measured  online  using  an  ultrasound  sensor  (Q401,  type 
LiquiSonic30  by  SensoTech,  GmbH,  Germany),  based  on  the 
influence  of  methanol  concentration  on  the  speed  of  sound  in 
methanol/water  solutions.  This  sensor  is  used  in  a  methanol 
concentration  controller,  which  uses  two  dosing  pumps  for 
pure  methanol  (P351)  and  pure  water  (P301)  (mzr-2905  by 
HNP  Mikrosysteme,  Germany)  as  actuators.  Methanol  con¬ 
centrations  between  zero  and  1.5  mol  dm-3  can  be  detected 
and  controlled.  The  flow  rates  of  the  dosing  pumps  can  be 
controlled  in  the  range  from  0.2  up  to  18  cm3  min-1.  To  ad¬ 
just  the  anode  pressure  and  also  to  strip  off  carbon  dioxide,  the 
recycle  vessels  are  equipped  with  a  nitrogen  purge/blanket. 
The  anode  pressure  can  be  controlled  in  the  range  between 
ambient  and  5  bars  absolute  (1.5  x  105  Pa).  The  liquid  inlet 
temperature  (which  is  also  the  DMFC  temperature  due  to  the 
applied  high  flow  rates)  can  be  controlled  in  the  range  be¬ 
tween  253  and  423  K.  The  DMFC  is  electrically  connected 
to  a  potentiostat  (HP60-50  by  Wenking  GmbH,  Germany), 
which  enables  operation  of  fuel  cells  from  below  1  W  up  to 
1  kW  at  a  maximum  of  50  A.  Galvanostatic  as  well  as  poten- 
tiostatic  operation  is  possible,  with  the  option  to  automati¬ 
cally  run  user-defined  load  scenarios. 

The  array  of  sensors  around  the  DMFC  enables  full  online 
material  balancing  of  all  key  components  (oxygen,  methanol, 
water,  carbon  dioxide). 


3.  Step  changes  in  cell  current 

3.1.  Experimental 

First,  a  steady  state  polarisation  curve  was  recorded  with 
the  following  operating  conditions: 

•  Anode: 

feed  temperature:  333  K  (equals  whole  DMFC  tempera¬ 
ture), 

methanol  feed  concentration:  1  mol  dm-3, 
flow  rate:  0.5  dm3  min-1, 
pressure:  1.7  bar; 

•  Cathode: 

feed  temperature:  293  K, 

feed  composition:  dry  air  (dew  point  ~27 1  K), 

flow  rate:  0.5  scbm  h- 1 , 

pressure:  1.7  bar. 

As  can  be  seen  in  Fig.  2,  the  OCV  lies  around  0.6  V,  and 
the  limiting  current  density  is  approximately  200  A  m-2.  The 
fairly  low  OCV  can  be  attributed  to  the  crossover  of  methanol 
and  is  in  accordance  with  literature  data  (0.65  V  in  [7],  0.55- 
0.7  V  in  [12]). 

At  temperatures  up  to  333  K,  it  is  known  that  methanol 
chemisorption  only  occurs  on  Pt  and  water  adsorption  on 


Fig.  2.  Recorded  polarisation  curve  with  illustration  of  conducted  current 
steps  (arrows)  and  the  activation-dominated  regime  (hatched  area,  left), 
pseudo-ohmic  regime  (middle)  and  transport-controlled  regime  (hatched 
area,  right). 

Ru,  respectively.  This  mechanism  changes  when  going  to 
higher  temperatures.  Then,  methanol  is  also  oxidised  on  Ru 
(see  [13]).  In  addition,  methanol  oxidation  accelerates  with 
higher  temperatures.  Since  the  aim  of  this  publication  is  the 
development  and  investigation  of  dynamic  DMFC  models, 
defined  conditions  and  reaction  mechanisms  for  modelling 
are  mandatory  for  obtaining  reasonable  results.  Hence,  an 
operating  temperature  of  333  K  was  chosen.  The  relatively 
low  performance  of  the  presented  DMFC  mainly  results  from 
this  moderate  operating  temperature. 

In  the  polarisation  curve  (Fig.  2),  the  activation-controlled 
and  mass-transport  controlled  regimes  are  marked  by  hatch¬ 
ing.  In  between,  without  hatching,  is  the  pseudo-ohmic 
regime.  When  investigating  cell  responses  to  current  steps, 
the  initial  current  as  well  as  the  final  current  should  strongly 
influence  the  voltage  response.  Hence,  a  series  of  current 
steps  between  the  different  regimes  were  performed.  These 
steps  are  visualised  by  arrows  in  Fig.  2.  Figs.  3-5  show  volt¬ 
age  responses  to  the  different  current  steps.  Each  of  the  cur¬ 
rent  densities  was  maintained  for  10  min  before  it  was  again 
set  to  the  previous  current  level.  This  assures  that  the  cell 
reaches  a  steady  state  before  the  next  current  change  is  ap¬ 
plied.  Voltage  responses  of  systems  which  are  not  in  steady 
state  prior  to  the  current  step  do  not  have  defined  conditions, 
and  interpretation  and  modelling  would  be  difficult. 

Fig.  3  shows  voltage  responses  to  steps  below  520  A  m-2. 
The  solid  line  shows  cell  voltage  responses  to  current  density 
steps  between  520  and  260  Am-2.  The  voltage  response  is 
quite  symmetrical,  showing  instantaneous  relative  overshoots 
of  approximately  40%  whenever  the  current  is  switched  be¬ 
tween  the  levels,  i.e.  the  absolute  value  of  the  voltage  change 
directly  after  the  current  step  exceeds  the  absolute  steady  state 
voltage  difference  by  40%.  The  steady  state  is  reached  after  a 
certain  relaxation  time.  In  contrast  to  this,  voltage  responses 
to  current  density  steps  between  OAm-2,  i.e.  open  circuit 
voltage  (OCV),  and  520 Am-2  are  non-symmetrical.  The 
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Fig.  3.  Cell  voltage  as  a  function  of  time  for  current  steps  from  the  lower 
pseudo-ohmic  regime:  Current  steps  between  0  and  520  A  m  2  (dotted), 
steps  between  260  and  520  A  m~2  (solid). 

voltage  response  at  OCV,  i.e.  when  the  current  is  switched 
off,  shows  an  overshoot  of  approximately  80%,  and  a  relax¬ 
ation  time  of  up  to  250  s,  in  comparison  to  roughly  50  s  in  the 
formerly  discussed  current  steps.  Additionally,  the  shape  of 
the  OCV  voltage  response  shows  rather  S-shape  than  expo¬ 
nential  behaviour  (like  in  the  former  step).  At  a  current  den¬ 
sity  of  260  A  m-2  the  activation-controlled  regime  is  reached 
and  the  forward  reaction,  i.e.  the  oxidation  of  methanol  takes 
place.  At  OCV  the  overall  reaction  stops  immediately,  and 
an  equilibrium  which  is  not  clearly  defined  is  reached.  This 
equilibrium  depends  on  the  oxidation  state  of  the  catalyst  as 
well  as  on  the  equilibria  between  the  various  methanol  reac¬ 
tion  intermediates.  This  process  takes  much  more  time  than 
reaching  a  steady  state  at  electric  current. 


Fig.  4.  Cell  voltage  as  a  function  of  time  for  current  steps  from  the  up¬ 
per  pseudo-ohmic  regime:  Current  steps  between  0  and  1050  A  m~2  (dot¬ 
ted),  steps  between  260  and  1050  A  m-2  (black),  steps  between  520  and 
1050 Am-2  (grey). 


Fig.  5.  Cell  voltage  as  a  function  of  time  for  current  steps  from  the  mass- 
transport  controlled  regime:  Current  steps  between  0  and  1570  Am  2  (dot¬ 
ted),  steps  between  260  and  1570  A  m-2  (black,  solid),  steps  between  520 
and  1 570  A  m  2  (light  grey),  steps  between  1050  and  1 570  A  m  ~2  (dark 
grey). 

Similar  effects  can  be  observed  for  steps  from  higher 
current  density  levels,  see  Figs.  4  and  5.  Steps  to  OCV  result 
in  much  higher  overshooting  and  longer  relaxation  times.  In 
contrast  to  Fig.  3,  voltage  responses  to  a  current  step-up  show 
small  fluctuations.  Furthermore,  the  voltage  overshoots  and 
subsequently  decays  to  new  quasi-steady  state  levels.  These 
new  quasi-steady  state  voltages  may  be  as  low  as  the  over¬ 
shoot  voltages  obtained  directly  after  the  current  step.  Decay 
and  fluctuation  should  be  assigned  to  diffusion  limitation 
effects.  While  the  upper  current  density  in  Fig.  4  is  always 
1050  A  m-2,  which  is  close  to  the  mass  transport  dominated 
regime  but  not  within,  the  upper  current  density  in  Fig.  5  is 
well  within  this  region.  This  leads  to  the  following:  While  the 
voltage  response  to  the  current  step-up  in  Fig.  4  is  decaying 
to  a  steady  state  cell  voltage  near  the  overshoot  cell  voltage 
due  to  mass  transport  limitation  effects.  Fig.  5  shows  intense 
fluctuations  of  the  voltage  signal,  attributed  additionally  to 
the  onset  of  CO2  bubble  formation.  Furthermore,  it  can  be  ob¬ 
served  that  the  system  takes  in  the  order  of  minutes  to  reach  a 
new  steady  state,  if  at  all  in  the  given  time  interval.  Apart  from 
these  effects,  when  stepping  down  to  lower  currents,  Figs. 
4  and  5  show  smooth  voltage  responses  with  exponential 
shapes  after  the  overshoots.  This  behaviour  is  comparable 
to  the  current  density  steps  from  520  to  260 Am-2  in 
Fig.  3. 

Fig.  6  plots  the  absolute  overshoots  of  cell  voltage  re¬ 
sponses  vs.  the  initial  and  final  cell  current  density  for  the 
conducted  current  steps.  Only  stable  voltage  signals  are  eval¬ 
uated.  As  can  be  seen,  increasing  final  cell  current  density 
correlates  with  decreasing  absolute  overshoot  (for  constant 
initial  cell  current  density).  An  already  mentioned  reason 
for  this  kind  of  behaviour  is  the  decay  of  cell  voltage  to  a 
lower  steady  state  voltage  after  an  overshoot,  which  can  be 
attributed  to  mass  transport  phenomena. 
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Fig.  6.  Voltage  response  to  cell  current  steps:  Absolute  overshoot  vs.  initial 
and  final  cell  current  density  (only  stable  voltage  signals  were  evaluated, 
dashed  lines  separate  the  regimes  in  the  polarisation  curve,  compare  with 
Fig.  2). 

Furthermore,  except  for  current  density  steps  from 
0  A  m-2  to  higher  current  densities,  higher  initial  cell  current 
densities  result  in  higher  overshoots  (for  constant  final  cell 
current  density).  This  leads  to  the  following  effect:  the  final 
steady  state  cell  voltage  is  mostly  independent  from  the  ini¬ 
tial  cell  current  density.  Hence,  the  cell  voltage,  the  methanol 
concentration  on  anode  and  cathode  side,  and  the  methanol 
crossover  should  reach  the  same  steady  state  after  a  current 
step  to  the  same  final  cell  current  density.  In  contrast  to  this, 
the  initial  conditions  vary  strongly.  Directly  after  the  cur¬ 
rent  step,  the  cell  voltage  is  strongly  influenced  by  the  initial 
steady  state  and  methanol  crossover.  Methanol  crossover  will 
be  lower  with  increasing  initial  cell  current  density.  This  may 
cause  lower  cathodic  overpotentials,  and  as  a  result  a  higher 
cell  voltage,  when  switching  to  the  same  final  cell  current 
density  from  higher  initial  cell  current  densities. 

Since  the  methanol  crossover  is  proportional  to  the  CO2 
concentration  on  the  cathode  side,  there  is  a  dependence 
of  the  CO2  concentration  on  current  density.  Schultz  et  al. 
[14,10]  have  shown  this  in  steady  state  measurements  and 
simulations.  Due  to  online  detection  of  CO2  at  the  cathode 
outlet  (after  the  condensor,  see  Fig.  1),  the  experimental 
setup  enables  dynamic  recording  of  the  CO2  concentration. 
Fig.  7  presents  the  CO2  content  in  the  cathode  exhaust  air 
as  a  function  of  time  for  the  current  steps  between  260  and 
1050  A  m-2.  Although  the  response  cannot  be  quantitatively 
interpreted  due  to  the  residence  time  influence  of  the  fixtures 
(pipes,  condensor).  Fig.  7  shows  qualitatively  the  response 
of  methanol  crossover  to  the  current  steps:  After  a  delay 
time  of  ~30  s,  the  CO2  concentration  gradually  moves  to  a 
higher  steady  state  when  stepping  down  to  260  A  m  ,  or  to 
a  lower  steady  state  when  stepping  up  to  1050  Am-2.  Here, 
a  higher  current  density  causes  higher  anodic  methanol 
conversion.  This  lowers  the  anodic  methanol  concentration, 
hence  the  diffusive  methanol  crossover  and  the  cathodic 
CO2  production.  The  opposite  holds  for  a  decrease  in  current 
density. 


Fig.  7.  Cell  voltage  (solid)  and  CO2  content  in  cathode  exhaust  air  (dashed) 
as  a  function  of  time  for  current  steps  between  260  and  1 050  Am  7 

3.2.  System  theoretical  analysis  of  process  models  of 
different  complexity 

3.2.1.  Set  of  governing  equations 

To  describe  and  analyse  the  behaviour  of  the  DMFC  with 
respect  to  dynamic  changes  in  current  density,  a  non-linear 
mathematical  model  was  developed  accounting  for  the  fol¬ 
lowing  phenomena  (see  Fig.  8): 

•  Diffusive  mass  transport  through  the  anode  diffusion  layer. 

•  Electrochemical  oxidation  of  methanol  in  the  anode  cata¬ 
lyst  layer. 

•  Formation  of  an  adsorbed  intermediate  COx  during  elec¬ 
trochemical  methanol  oxidation. 

•  Electrochemical  reduction  of  oxygen  in  the  cathode  cata¬ 
lyst  layer. 

•  Methanol  crossover,  i.e.  undesired  transfer  of  methanol 
through  the  PEM  by  diffusion  and  electroosmosis. 

•  Undesired  electrochemical  oxidation  of  methanol  at  the 
cathode  catalyst  layer. 

The  following  model  assumptions  were  made: 

•  Residence  time  in  the  anode  compartment  is  very  low,  so 
that  the  compartment  effectively  is  a  differential  reactor. 
Hence,  the  concentration  in  the  compartment  is  assumed 
to  be  equal  to  the  inlet  concentration. 

•  Ohmic  drops  in  current  collectors  and  electric  connections 
are  negligible. 

•  The  fuel  cell  is  operated  isothermally. 

•  Oxygen  is  fed  in  excess,  i.e.  oxygen  conversion  in  the  cath¬ 
ode  compartment  is  negligible,  and  therefore  no  oxygen 
mass  balance  is  required. 

•  Oxygen  and  carbon  dioxide  do  not  diffuse  into  the  PEM. 

•  Water  concentration  is  assumed  to  be  constant  (excess 
component  in  a  liquid  mixture). 

•  Mass  transport  resistances  in  the  catalyst  layers  are  negli¬ 
gible  due  to  the  fact  that  these  are  thin  (35  pun)  in  compari¬ 
son  to  the  diffusion  layer  (170  pom)  and  the  PEM  (100  pun). 
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Fig.  8.  Illustration  of  the  DMFC  model  including  reactants  and  physico¬ 
chemical  phenomena  (A:  anode  compartment;  AD:  anode  diffusion  layer; 
AC:  anode  catalyst  layer;  M:  membrane;  CC:  cathode  catalyst  layer). 


lumped  intermediate  COx  is  introduced.  While  the  first  re¬ 
action  step  is  the  reaction  of  methanol  to  the  intermediate, 
the  second  reaction  is  the  oxidation  of  this  intermediate  plus 
an  integrated  reaction,  which  is  the  adsorption  of  water  on 
Ru.  This  adsorption  step  on  Ru  is  known  to  be  fast  compared 
to  the  slow  methanol  oxidation,  see  [20],  and  therefore  it  is 
assumed  to  be  in  equilibrium. 

As  can  be  seen  from  the  assumptions,  the  model  is  lim¬ 
ited  to  an  operating  regime  of  very  high  anodic  volume  flow, 
a  lumped  reaction  mechanism  and  the  lower  pseudo-ohmic 
regime  of  the  polarisation  curve.  Under  these  conditions,  an 
analysis  of  the  dynamic  processes  occurring  inside  the  MEA 
can  be  conducted  with  minimal  influence  of  mass  transport. 

Based  on  the  assumptions  given  above,  balance  equations 
for  methanol  concentration  in  the  catalyst  layer  (Eq.  (3)), 
COx  surface  coverage  (Eq.  (4)),  methanol  concentration  in 
the  membrane  (Eq.  (6)),  anodic  overpotential  (Eq.  (5))  and 
cathodic  overpotential  (Eq.  (7)),  as  well  as  reaction  rates  (Eqs. 
(8)-(10))  were  formulated: 


•  Mass  transport  in  the  diffusion  layers  is  quasi-steady  state. 

•  On  the  anode  side,  a  pure  liquid  phase  mixture  is  assumed, 
i.e.  gas-phase  formation  by  release  of  carbon  dioxide  bub¬ 
bles  is  not  taken  into  account.  Problems  with  bubble  for¬ 
mation  only  get  visible  when  working  and  modelling  in 
or  near  the  transport-controlled  regime  of  the  polarisation 
curve,  which  is  not  done  here. 

•  Methanol  which  is  transported  to  the  cathode  is  assumed 
to  be  instantaneously  reacting  with  oxygen  at  the  cathode, 
hence  its  concentration  in  the  cathode  catalyst  layer  is  zero. 

•  Experiments  were  carried  out  with  zero  pressure  differ¬ 
ence  between  anode  and  cathode  compartment,  hence 
no  pressure-driven  transport  through  the  membrane  takes 
place. 

•  In  the  modelled  scenarios,  all  reactions  are  far  from  equi¬ 
librium.  Hence,  Tafel  kinetics  are  used  for  description. 

•  Due  to  the  postulated  instantaneous  dissolution  of  carbon 
dioxide  after  production  and  due  to  the  applied  Tafel  ki¬ 
netics,  no  carbon  dioxide  balance  is  formulated. 

•  A  two  step  reaction  mechanism  is  used  to  describe  the 
electrochemical  oxidation  of  methanol  at  the  platinum- 
ruthenium  catalyst.  The  mechanism  is  a  lumped  approach 
and  accounts  for  the  oxidation  of  a  Pt-adsorbed  interme¬ 
diate  COx  (e.g.  CO): 

CH3OH  +  Pt  — ►  Pt-COx  +  4H+  +  4  e“  (1) 

Pt-C0x  +  H20^C02 +2H+ +  2e“ +Pt  (2) 

There  is  intense  discussion  about  the  real  methanol  oxi¬ 
dation  mechanism  and  the  intermediates  involved,  see  [15- 
19].  Focus  of  this  work  is  not  on  determining  or  validating 
a  specific  reaction  mechanism.  Instead  the  fact  is  used  that 
methanol  oxidation  is  known  to  occur  via  a  carbon  interme¬ 
diate  which  is  adsorbed  on  active  Pt  surface  sites  and  which 
tends  to  block  the  sites.  As  this  intermediate  may  be  CO, 
COH,  CHO  or  other  postulated  and  measured  species,  the 
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Transport  of  methanol  through  the  membrane  is  described 
by  Eq.  (1 1),  consisting  of  a  Fickian  diffusion  term  (first  term 
on  r.h.s.),  and  an  electroosmotic  term  (second  term)  with  the 
lumped  electroosmotic  parameter  k^2.  The  latter  is  derived 
by  using  the  isobaric  convective  flow  velocity  equation  given 
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by  Schloegl  [21]: 

„M  _  _rM  r#  m\ 

^CH3 OH, convective  cCH30H  ^  ^  v1-5/ 

Analogous  to  Sundmacher  et  al.  in  [7],  the  electrostatic  gra¬ 
dient  dcp/dz  can  be  formulated  in  terms  of  the  cell  current 
density  z'ceii,  using  the  Nernst-Planck  equation  for  proton 
transport: 


_  (cell 

dF  =  F2(D“  /RT  +  c“  k+hi) 


(14) 


which  is  inserted  into  Eq.  (13).  From  this,  a  lumped  elec- 
troosmotic  parameter  is  obtained: 


^02 


k(j) 

/RT  +  c&  4//x) 


05) 


The  overall  cell  voltage  C/Ceti  (Eq.  (12))  is  a  function  of  the 
standard  open-circuit  cell  voltage  J/q,  the  overpotentials  at 
anode  and  cathode,  rj a  and  ;?c,  respectively,  and  the  Ohmic 
loss  in  the  PEM,  Re\  =  dM /km.  The  overpotentials  are  de¬ 
fined  with  regard  to  the  standard  electrode  potentials. 


3.2.2.  Linear  system  analysis  and  overshoot 

After  linearisation,  the  given  set  of  balance  equations.  (3)- 
(7),  can  be  cast  into  the  following  standardised  state  space 


form: 

dx 

—  =  Ax(t )  +  Bu(t ) 
df 

(16) 

y(t)  —  Cx(t)  +  Du(t) 

(17) 

where  A,  B,  C  and  D  are  the  system  matrix,  input  matrix, 
output  matrix  and  direct  transmission  matrix,  respectively,  u 
the  perturbed  input  variable,  i.e.  in  this  case  the  cell  current 
density,  y  the  observed  output  variable,  i.e.  the  cell  voltage 
response,  and  x  is  the  vector  of  state  variables. 

Laplace  transformation  (t  — »■  s;  y(t)  — ►  Y(s)\  u(t)  -> 
U(s)\x(t)  — »■  X(.v))  of  this  equation  set  yields  the  transfer 
function  G,  which  is  the  ratio  of  output  to  input  variable: 

G(s)=ffr  =  B(sI-A)-lC  +  D  (18) 

U(s) 

Y (s)  is  the  Laplace  transformed  cell  voltage,  and  U(s)  is  the 
Laplace  transformed  cell  current  density.  Instead  of  calculat¬ 
ing  the  global  transfer  function  using  all  equations  at  once, 
the  system  can  be  decomposed  into  single  modules  which 
are  interconnected  as  shown  in  the  blockdiagram  in  Fig.  9. 
This  decomposition  is  done  with  regard  to  the  different  func¬ 
tional  layers  of  the  DMFC:  Anode  catalyst  layer  (G a),  cath¬ 
ode  catalyst  layer  (Gc)  and  membrane  (Gm).  Subsequently, 
the  anode  catalyst  layer  is  decomposed  into  the  different  state 
variables:  Methanol  concentration  (G acX  COx  surface  cov¬ 
erage  (Gat)  and  anode  catalyst  charge  (Gae)-  This  method 
of  decomposition  has  the  advantage  that  the  mathematical  de¬ 
scription  of  the  single  modules  can  be  chosen  from  within  a 
library.  Besides  the  descriptions  derived  from  Eqs.  (3)— ( 12), 
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Fig.  9.  Representation  of  the  DMFC  as  a  transfer  function  block  diagram 
(left)  and  an  overview  of  the  presented  model  variants  (right):  G  =  const 
quasi-steady  state  assumed. 


this  library  contains  various  alternative  mathematical  descrip¬ 
tions  for  each  state  variable.  The  appropriate  description  has 
to  be  chosen  according  to  the  model  variant.  Four  different 
model  variants  are  investigated  in  the  following  chapters,  an 
overview  is  given  in  Fig.  9  (r.h.s.). 

The  experimental  input  u(t )  shows  a  step  in  cell  current 
density,  which  in  the  Laplace  domain  gives  U(s)  =  Au/s , 
where  Am  is  the  step  size  Am  =  zCcii.r>o  —  z’ceii.ss-  The  index 
ss  always  refers  to  the  steady  state  prior  to  the  current  step. 

Based  on  the  observation  that  the  experimental  voltage  re¬ 
sponses  typically  overshoot  directly  after  the  current  step,  a 
criterion  in  the  Laplace  domain  for  overshooting  is  formu¬ 
lated:  The  transfer  function  G(s)  is  written  in  the  pole-zero 
form: 


G(s)  =  k 


n:u(* 


n:u(* 


Z/) 

Pi) 


(19) 


where  Z,  are  the  zeros  or  roots  of  the  numerator  polynomial 
and  Pi  are  the  poles,  i.e.  the  roots  of  the  denominator  polyno¬ 
mial.  The  number  of  poles  and  zeros  is  the  same,  as  soon  as 
D  ^  0.  As  illustrated  in  Fig.  10,  in  case  of  overshooting,  the 
system  response  has  its  maximum  absolute  value  at  t  — »•  0+, 
i.e.  lim^o  |y(f)l<  and  decreases  for  t  ->  oo  towards  a  new 
steady  state  with  a  smaller  absolute  value,  i.e.  lim^oo  |y(f)|. 
Final  value  theorem  and  initial  value  theorem  then  result  in 


lim  y(t)  =  lim  T(y)s  =  lim  G(s)Au  =  kAu  (20) 

t—>  0  S—>  oo  s—>  oo 


lim  y(t)  —  lim  T(j)s  = 

r->oo  5— >0 


lim  G(s)Au  =  kAu 

s— >o 


n;-Li(-z,) 

n;=i(-^o 


(21) 


And  the  condition: 


lim  |  v(0 1  >  lim  |y(r)| 

r->  0  t~ro o 


(22) 
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Fig.  10.  Principle  of  overshooting  systems. 


gives  the  overshoot  criterion: 

n  n 

ii  n  ^ii  >  ii  n  z/ii  (23) 

i=i  i=t 

The  following  chapters  investigate  four  model  variants. 


Accordingly,  the  number  of  exchanged  electrons  in  Eq. 
(5)  has  to  be  changed  from  4  to  6,  and  Eq.  (4)  is  dismissed. 

Also,  the  description  of  mass  transport  through  the  mem¬ 
brane  is  simplified.  The  most  simple  mass  transport  model  is 
a  Fick  diffusion  approach,  assuming  a  negligible  electroos- 
motic  drag,  hence  Eq.  (11)  reduces  to 


Dm 

„M  _  UC\hO\\  f\C 
nCH3OH  —  CCH3OH 


(26) 


This  means  that  a  linear  concentration  profile  in  the  mem¬ 
brane  is  assumed,  and  thus  the  methanol  flux  density  of 
methanol  is  constant  all  over  the  membrane: 


M.A  _  M.C  _  M 
nCH3OH  —  nCH3OH  —  "ClhOII 


(27) 


3.2.3.  Reaction-model  with  one  state  variable 
( lx-model) 

The  given  set  of  balance  equations  (3) — (7)  represents  dif¬ 
ferent  dynamic  processes  which  interact  with  each  other.  To 


Fig.  9  (r.h.s.)  gives  an  overview  of  the  transfer  functions 
for  the  lx-model.  Calculation  yields  the  following  trans¬ 
fer  functions,  where  in-  and  outputs  are  in  the  same  or¬ 
der,  from  top  to  bottom,  as  in  the  blockdiagram  in  Fig.  9 
(l.h.s.): 


G  c 


study  and  understand  the  model  and  the  influences  of  the 
various  phenomena,  a  study  of  the  individual  dynamic  pro¬ 
cesses  is  helpful.  A  minimum  dynamic  model  with  one  state 


Gm  = 


^  0  'Dch3oh/^M 

dM/KM  0 
\  0 


/ 


-As/Vac 


GacG)  = 


S  +  £AD  As/  VAC  +  As/  VACrAio  exp(aA1  Fr,AiSS/R/  T)/c^f0H 
~  As/  VACr  AioaAi  e/  R/  T  exP  (aAi  FrjAiSS/ R/  7")cch3oh,ss/ cch3oh 
V  *  +  *A%/  ^Ac  +  Ag/  VACrA1o  exp(«A,  Ftja^/R/  T)/c^OH  ) 


(28) 


(29) 


(30) 


variable,  denoted  in  the  following  as  lx-model,  is  derived 
from  the  set  of  Eqs.  (3)-(12).  The  most  influential  state  vari¬ 
able  is  the  methanol  concentration  in  the  catalyst  layer:  It  is 
changing  during  a  current  step,  influencing  anode  and  cathode 
overpotential,  as  well  as  the  methanol  profile  within  the  mem¬ 
brane.  For  the  lx-model  the  following  assumptions  are  made: 


•  The  balance  equations  for  charge,  Eqs.  (5)  and  (7),  and 
methanol  inside  the  membrane,  Eq.  (6),  are  in  a  quasi¬ 
steady  state. 

•  The  balance  equation  for  methanol  inside  the  anode  cata¬ 
lyst  layer,  Eq.  (3),  is  the  only  dynamic  equation. 

•  The  methanol  oxidation  is  assumed  to  take  place  without 
the  reaction  intermediate  COx,  therefore  the  rate  expres¬ 
sion  rAi  has  to  be  modified  to  match  the  overall  reaction: 

CH3OH  +  H20  C02  +  6H+  +  6  e“  (24) 


and 


..AC 

tCH3OH 

r\\  =  CAIO-J^C - exp 

cCH3OH 


ot\\F  ) 
RT 


(25) 


GAe(a)  = 


RT 


RT 


“A1  GcA^3oH  , 


«A1  Cf'celt' 


(31) 


Applying  the  overshoot  criterion,  Eq.  (23),  on  the  derived 
model  equations  results  in  the  following  overshoot  criterion: 


6  FI), 


tmax  — 


'M 
CH3OHa 


;.AD  _AC 


“CH3OH.ss 


(aAi  +ac  -  1) 


rfM(/tAD(l  -  ac)  +  kMaAl) 


(32) 


Obviously,  there  is  a  maximum  initial  steady  state  current 
density  imax,  above  which  in  any  case  no  overshooting,  and 
below  which  always  overshooting  will  occur.  It  should  be 
emphasized  that  with  such  a  simple  model  of  only  one  state 
variable,  the  voltage  response  can  already  show  overshoots 
depending  on  the  parameters.  In  case  that  both  charge  trans¬ 
fer  coefficients,  aA\  and  ac,  are  0.5,  ;max  is  0,  therefore  no 
overshoots  can  occur.  Fig.  1 1  (top)  shows  the  experimental 
cell  voltage  response  for  a  current  density  step  from  520  to 
260  Am-2.  In  order  to  simulate  this,  all  known  parameters 
(values  see  Nomenclature)  are  inserted  in  the  model  and  the 
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charge  transfer  coefficients  «ai  and  ac  are  set  to  0.5  (stan¬ 
dard  value).  Parameters  unknown  are  the  anodic  and  cathodic 
reaction  rate  constants,  r,\  i  o  and  rco,  respectively.  They  are 
estimated  according  to  the  following  criteria.  First,  the  ca¬ 
thodic  reaction  should  be  much  faster  than  the  anodic  reac¬ 
tion.  Second,  the  experimentally  observed  polarisation  curve 
(Fig.  2)  has  to  be  fitted.  Accordingly,  the  cathodic  reaction 
rate  constant  is  chosen  to  be  two  orders  higher  than  the  an¬ 
odic  reaction  rate  constant,  and  the  reaction  parameters  are 
adjusted  to  fit  the  experimental  polarisation  curve.  Fig.  1 1 
(bottom,  dotted  line)  shows  the  simulated  voltage  response 
(reaction  parameters  are  given  ibidem).  With  this  set  of  pa¬ 
rameters,  no  overshooting  can  occur  according  to  Eq.  (32). 
The  voltage  response  for  charge  transfer  coefficients  >0.5 
is  shown  in  Fig.  1 1  (bottom,  solid  line).  As  expected  for  the 
given  parameter  set,  the  voltage  response  shows  overshooting 
Oceli.ss  <  'max)-  Qualitatively,  the  model’s  voltage  response 
is  similar  to  the  experimental  voltage  response.  Both  have  a 
maximum  at  t  =  0  s  and  decrease  slowly  to  the  new  steady 
state.  Nonetheless,  quantitative  comparison  shows  significant 


ometry  and  the  layers’  morphology,  other  additional  storages 
have  to  exist. 

3.2.4.  Reaction-adsorption  model  with  two  state 
variables  (2x-model) 

As  the  1  A-model  cannot  describe  the  experimental  obser¬ 
vations  appropriately,  the  model  is  extended  by  introducing 
a  second  state  variable  Ocox ,  hence  another  storage.  As  ex¬ 
plained  already  in  Section  3.2.1,  it  is  known  that  the  elec¬ 
trochemical  methanol  oxidation  follows  a  multi-step  mecha¬ 
nism  including  at  least  one  adsorption  step,  e.g.  like  in  Eqs. 
(1)  and  (2).  The  mass  balance  of  the  intermediate  step  is 
given  in  Eq.  (4),  with  the  reaction  rates  r^i  and  r,\2  as  in 
Eqs.  (8)  and  (9).  The  description  of  mass  transfer  through  the 
membrane  is  identical  to  the  one  in  the  1  A-model.  The  new 
model  will  be  referred  to  as  2x-model  due  to  its  two  state 
variables. 

As  in  the  I  A-model,  the  equation  system  is  transformed 
into  Laplace  domain.  G c  and  Gm  are  identical  to  the  ones  in 
the  lA-model.  Due  to  the  new  state  variable  dcox .  the  transfer 
functions  for  the  anode  (in-  and  outputs  are  in  the  same  order 
as  in  Fig.  9)  change  as  follows: 


-AS/VA 


s+kAOAs/VAC+As/VAcrA10txp(aAFr1Kss/R/T)/c^OHa-ecOx.ss)/(l-e^Ox) 

x) 


As/  VAcrA10  exp(aAFijA,ss/.K/7>cH  ,OH  ,ss/c! 


3 
.ref 
CH3OH 


A1- O 


GAC(S)  -  i+*ADAs/yAC+As/yACrAloexp(aA^Ass/S/r)/(^f3OH(1_0COxss)/(1_e«fx) 

-As/VAcrA10aAF/ff/rexp(aAF^AjSS/ff/r)cA|3OHiSS/^;iOH(l-ecox.ss)/(l-e^x) 
s+kADAs/VAC+As/VAcrAioeiLp(aAF^ss/R/T)/(^1  0H(l-<fcOx,ss)/(l-@cox)  J 


Gat(s) 


/  aAF/R/T exp(o'AFf)AiSS/i?/r)/nmax,pt(/-AIOfAg3OH  ss/c^3OH(l-ecOx,ss)/(l-tlcox)_''A2O0COx,ss/tlceoxl  \ 

s-l/nmMM-rAi0exp(aAFriA'SJR/T)c^OH'SJc,l!*l3O}1/G-8cbx,s!)-rA20exP(o,AFvA,ss/R/T')/9,l!;f3x) 
rAioexp(aAF)jAiSS/it/r)/cgl3OHa-ecOx,ss)/(l-0cox)/',max,Pt 
V  t /"max.Ptt- '  A10  Oxpfo-A  GA,xs/'t/'r)<'Q[.A)j[  ^s/rf  F[,on  A  ^  “^COx,ss  )-rA20  exp (otAFi]A,ss/ R/ D/^cOx)  ) 


Gae(s)  = 


/ _ -jtr(-2>Aio9gfx+2rAioegfxeCox.ss) _ 

(2'A10ccH3OH,ss®COx^c°x.ss—  1  )+'W2O0COx,sscgH3oH^COx,ss_  B)“A^ 

a>j?r(2rAi0cg§3OHssegQx— rA20cgH3OH+''A20cCH3OH6lCOx) 

(2<  Alocg{j30H  ss0gox(0COx,ss-l)+''A2O0COx,ssegg3QH(0gQx— lfio-Af 
RT 

01 A  ^cell.ss 


(33) 


(34) 


(35) 


differences.  The  model’s  absolute  overshoot  and  the  relax¬ 
ation  time  are  much  smaller  than  the  experimental  values. 
The  absolute  overshoot,  i.e.  the  difference  between  maximum 
voltage  and  new  steady  state  voltage,  is  20  mV  in  the  experi¬ 
ment  and  approximately  10  times  smaller  in  simulation.  This 
may  originate  from  the  selected  set  of  reaction  parameters. 
The  more  important  deviation  is  the  relaxation  time:  while 
in  the  experiment  it  takes  approximately  40-50  s  to  reach  a 
new  steady  state,  the  steady  state  is  reached  within  5  s  in  the 
model.  This  indicates  that  the  chosen  storage  capacity  within 
the  model  is  insufficient.  As  the  only  storage  in  the  model  is 
the  anode  catalyst  layer  volume  which  is  given  by  cell  ge- 


where  the  anodic  charge  transfer  coefficients  are  assumed 
to  be  identical  (a a  =  aAi  =  o;a2)-  The  overshoot  criterion 
can  be  derived  analogously  to  the  l.A-model.  All  previous 
parameters  are  used,  except  for  the  reaction  rate  constants 
aaio  and  r,\20-  The  cathodic  reaction  rate  constant  of  the 
1  A-model,  including  the  assumption  of  a  faster  cathodic 
reaction,  is  adopted.  The  criterion  for  overshoot  yields 
an  implicit  equation  with  three  variables:  r,\  i o •  ''A20  and 
the  initial  steady  state  current  density  fCell,ss-  Since  this 
equation  is  lengthy  and  complicated,  it  is  not  presented 
here.  Furthermore,  it  cannot  be  solved  explicitly,  but  can 
be  plotted  as  a  plane  in  M3  with  an  algebraic  mathematics 
program  like  Maple,  as  shown  in  Fig.  12  (top),  where  iceii,ss 
is  plotted  vs.  the  two  reaction  constants.  The  plane  separates 
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Fig.  1 1 .  Experimental  voltage  response  (top)  and  voltage  response  of  the  lx- 
model  (bottom)  to  a  current  step  from  520  to  260  A  in  2 :  Overshoot  response 
with  i7a  i  =  ac  =  0.6,  r  a  i  o  =  10  mol  m  2  s_1  (solid),  non-overshoot  re¬ 
sponse  with  «ai  =  ac  =  0.5,  7a  i 0  =  1.4  x  10  mol  m  2  s_1  (dotted). 

the  overshoot  region  (left)  from  the  non-overshoot  region 
(right).  It  can  be  seen  that  the  plane  shows  an  incline.  At 
each  cell  current  density,  there  is  a  ratio  of  caio  versus  rA2o 
above  which  overshooting  occurs.  This  ratio  is  increasing 
with  increasing  cell  current  density.  Hence,  similar  as  in  the 
lx- mode  I,  for  a  set  of  reaction  rate  constants  a  maximum 
current  density  exists,  above  which  no  overshoot  can  occur. 

Such  behaviour  was  not  observed  in  the  experiments  (see 
Fig.  5),  steps  from  high  initial  cell  current  densities  still  pro¬ 
duce  overshoots.  Mass  transport  limitation,  which  in  the  mod¬ 
els  is  only  accounted  for  in  a  simplified  form,  in  the  experi¬ 
ments  leads  to  a  lower  limiting  current  density.  This  current 
density  may  be  lower  than  the  model’s  maximum  current  den¬ 
sity  for  overshooting,  (’cell, max,  hence  according  to  the  model 
overshoot  will  always  occur.  Furthermore,  the  cell  current 
steps  conducted  at  higher  current  densities  are  far  from  be¬ 
ing  linear  system  deflexions  (see  Section  3.2.6),  for  which 
exclusively  the  model  holds. 

Fig.  12  (bottom)  shows  the  separation  line  for  the  experi¬ 
mental  initial  current  density  of  520  A  m-2.  Here,  sets  of  pa¬ 
rameters  in  the  lower  right  region  will  result  in  overshooting 
behaviour.  Due  to  the  experimentally  observed  overshooting, 
the  anodic  reaction  parameters  should  be  within  the  over¬ 
shooting  regime  of  Fig.  12  (bottom).  The  anodic  reaction 
parameters  are  fitted  to  the  experimental  polarisation  curve, 
and  the  dynamic  answer  to  the  current  density  step  from  520 
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Fig.  12.  Analytically  calculated  separation  plane/line  between  overshooting 
and  non-overshooting  region:  Anodic  reaction  constants  vs.  initial  current 
density  (top),  anodic  reaction  constants  at  initial  current  density  (Ceii,ss  = 
520  A  m~2  (bottom). 

to  260  A  m~2  is  simulated.  The  result  is  shown  in  Fig.  13 
(parameters  are  given  ibidem  or  in  Nomenclature).  While  the 
absolute  overshoot  is  still  ten  times  lower  than  the  experimen¬ 
tal  one,  the  relaxation  time,  approximately  40  s,  is  close  to  the 
experimental  relaxation  time,  indicating  more  realistic  mod¬ 
elling  results.  A  different  selection  of  reaction  parameters 
may  lead  to  higher  or  lower  overshoot,  but  will  not  influence 
the  relaxation  time.  Therefore,  a  significant  step  towards  a 
good  fit  to  the  experiments  is  achieved  by  the  2x-model. 


Fig.  13.  Voltage  response  of  the  2x-model  to  a  current  step  from  520  to 
260  Am-2  with  taio  =  1 0-6  mol m~2  s_  1 ,  r^io  =  5  x  10-8  molm-2  s-1. 
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Fig.  14.  Concentration  and  overpotential  responses  of  the  non-linear  2x- 
model  to  a  current  step  from  520  to  260  A  m-2.  From  top  to  bottom:  COx- 
coverage,  methanol  concentration  in  the  anode  catalyst  layer,  anodic  over¬ 
potential.  cathodic  overpotential. 


Fig.  14  shows  the  result  of  the  simulation  of  the  numer¬ 
ically  evaluated  non-linear  2x-model  (see  Section  3.2.6).  It 
shows  the  transient  response  of  $cOx,  cch3oH’  3a  and  3C- 
As  can  be  seen,  the  anodic  overpotential  is  quite  insensi¬ 
tive  to  changing  concentrations  of  methanol  and  to  changing 
COx  surface  coverages  in  the  anodic  catalyst  layer  for  time 
t  >  0,  while  the  cathodic  overpotential  rjc  shows  a  signifi¬ 
cant  overshoot  which  is  the  reason  for  the  overall  cell  voltage 
overshooting.  The  theoretical  result  of  an  insensitivity  of  the 
anodic  overpotential  to  changes  in  methanol  concentration 
has  been  validated  in  half-cell  measurements  on  a  real  fuel 
cell  MEA  [22].  As  soon  as  the  methanol  concentration  in  the 
anode  catalyst  layer  reaches  a  steady  state,  the  methanol  flux 
through  the  membrane  reaches  a  steady  state,  and  with  this 
also  the  cathodic  overpotential. 

This  theoretical  result  supports  the  interpretation  of  the 
experimental  observations  at  the  end  of  Section  3. 1  and  Fig.  6. 

3.2.5.  Reaction-adsorption  model  with  two  state 
variables  and  a  distributed  membrane  (2.x  +  m-model) 

The  lx-model  and  the  2x-model  are  focussed  on  the  reac¬ 
tion  mechanism  and  its  influence  on  the  dynamic  behaviour 
of  the  DMFC.  The  membrane  dynamics  is  neglected  there 
by  assuming  a  steady-state  linear  concentration  profile  of 
methanol  in  the  membrane.  But  recent  studies  by  Schultz 
[10]  indicate  that  the  concentration  profile  of  methanol  in  the 
membrane  is  not  necessarily  linear,  and  the  membrane  is  a 
significant  storage  for  methanol.  To  investigate  the  influence 
of  the  dynamics  of  the  membrane  on  the  cell  voltage  response, 
the  2x-model  is  extended  with  a  further  state  variable,  which 
is  the  local  concentration  of  methanol  in  the  membrane.  This 
model  is  denoted  as  2x  +  m -model.  Compared  to  the  setup  of 
the  previous  transfer  function  elements  (G c,  G,\c (■'’),  G\t(.v) 
and  Gae  of  the  2x-model  hold  for  this  model,  too),  the  setup 


of  the  dynamic  membrane  transfer  function  element  Gm(s)  is 
more  complex,  since  mass  transport  through  the  membrane 
is  given  by  a  partial  differential  equation,  see  Eqs.  (6)  and 
(11).  The  procedure  is  as  follows.  As  can  be  seen  in  Fig.  9, 
the  input  into  the  membrane  module: 


GmG)  = 


Ks)  \ 
<-’ch3oh(^')  ) 


and  output  out  of  the  membrane 


ImC?)  = 


(  A$£oh (*)  \ 
dM/KM  ■  I(S ) 

xtM,A  /  \ 

\  ^C’lhOH^  / 


(36) 


(37) 


are  coupled  by  a  3  x  2  transfer  function  matrix  Gm(j): 


Yu(s)  =  GmW  •  UM(s) 

with 


/ Gmii(s)  GM12C5)' 


GmG)  = 


dM/icM  0 
\  Gm3i(s)  GM32G), 


(38) 


(39) 


GmiiC?),  Gmi2('5),  Gm3iCs)  and  Gm32(s)  are  functions  of  s 
and  are  determined  by  using  Eq.  (6),  which  is  linearised  at 
the  steady  state  operating  point. 

Assuming  the  following  two  conditions: 

•  methanol  concentration  in  the  membrane  at  the  spatial  co¬ 
ordinate  z  =  0,  where  the  membrane  is  in  contact  with  the 
anode  catalyst  layer,  is  identical  to  the  concentration  in  the 
anode  catalyst  layer,  i.e. 

cch3oh(-  =  0)  =  Ccn3OHi  (40) 

•  methanol  concentration  in  the  membrane  at  z  =  dM,  where 
the  membrane  is  in  contact  with  the  cathode  catalyst  layer, 
is  assumed  to  be  zero  due  to  the  total  oxidation  of  methanol 
with  oxygen,  i.e. 

cch3oh(£  =  ^M)  =  cch3oh  =  0;  (41) 

the  following  partial  differential  equation  results: 


dx 
-  =  D 


dt 


ch3oh 


dlx 


dx 


^2  -  42Gell,ss^  + 


i  2  AC 

*02  tCH3OH,ss‘cell,ss  e 


3  Z 

Pess(z/dM- 1) 


D, 


ch3oh 


1  —  e 


-Pe  sc 


Ocell  *cell,ss) 


(42) 


where  the  Peclet  number: 


Pe  ss 


^02^^cell,ss 

dm 

^CHsOH 


(43) 
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and  the  definition  a-  =  (cch3oh  —  Cch3oh  ss)  anc*  initial  con¬ 
dition  x(t  —  0,  z)  —  0  are  holding,  as  well  as  the  boundary 
condition: 

x(t,  z  =  0)  =  Cc§3oH  —  Cch3OH,ss’  x(?’  2  =  fl,M)  —  0 

(44) 


Laplace  transformation  yields  the  following  ordinary  dif¬ 
ferential  equation: 


Z)ch3ohx"('s>  z)  -  k^icen.ssS'is,  z )  -  sX(s,  z) 
AC 

cCH3OH,ssLell,ss  i 


^022CpH,nrr  «*ce  11, ss  ePe^(z/dM  1) 


Dm 
]  l3OI  f 


1 


*-Pes  s 


/(s) 


(45) 


with  X(j,  z  =  0)  =  C££30H(.s),  X(s,  z  —  dM)  =  0. 

After  solving  this  equation  analytically  and  inserting  the 
result  X(s,  z)  into  the  Laplace  transform  of  the  linearised  Eq. 
(11): 

N(s,  z)  =  -£>CH3OtA'(s,  z)  +  £02*cell,ssX(s,  Z) 

+  fc02CCH3OH,ss(z)!C?)  (46) 


a  Laplace  transform  of  the  methanol  mass  flow  is  obtained, 
which  can  be  cast  into  the  following  form: 

N(s,  z)  =  b\(s,  z)X(s ,  z)  +  b2(s,  z)I(s)  (47) 


b\  and  b2  are  complicated  functions  of  .v  and  z.  They  are  not 
presented  here  due  to  their  size. 

The  transfer  function  elements  demanded  in  Eq.  (39)  are 
obtained  from  Eq.  (47)  at  the  coordinates  z  =  0  and  z  —  dM: 

&022CcH3OH,ss*cell,ss  Q~Pe^!2  1 

2afM  1— e_Pess  s 

■  [~Pess  tPe^2  -  2 F&)  +  2ePe^2F2(s)]  (48) 


GmhW  = 


Gm12(s)  = 


dm 

^ch3oh 

dM 


ePess  F i  (s) 


(49) 


Unfortunately,  these  transfer  functions  are  transcendental 
functions,  and  using  them  in  this  form  would  not  result  in  a 
rational  transfer  function  for  the  fuel  cell  (Eq.  (19)),  hence 
neither  can  this  function  be  used  for  controlling  purposes 
nor  can  an  overshoot  analysis  be  carried  out  according 
to  Eq.  (23).  There  exist  several  methods  to  approximate 
transcendental  transfer  functions  by  rational  functions, 
see  e.g.  [23].  A  reliable  analytical  solution  is  obtained 
via  a  combination  of  the  residue  theorem  and  the  Laplace 
transformation.  Step  one  is  the  transformation  of  the  inverse 
Laplace  transformation  integral  into  a  contour  integral  with 
a  semicircle  of  radius  p  — »■  oo: 

1  rc+j°°  l  r 

y(t)  =  —  /  y (s)  e'7d.v  =  —  <£  Y(s)C‘ds  (54) 

J  Jc—jo O  Z7tJ  j  p-^-OQ 

If  all  singularities  of  Y(s)  are  inside  this  region  and  if 


lim  Y(s)  =  0  (55) 

s — >  OO 


according  to  the  residue  theorem  the  system  answer  in  the 
time  domain  y(t)  can  be  calculated  by  the  sum  of  the  contour 
integrals  around  each  singularity  P,: 


y(t)  = 


2jr;  Jp 


p— >oo 


Y(s)  ev,d.s  =  y  res,- 


i=  1 


(56) 


with 

res,-  =  — ;  i j)  T(s)esfdj-  (57) 

2 7tj  J o,  at  Pi 

Integration  of  this  equation  yields: 

res,-  =  lim  (.?  —  P,)Y(s)est  =  const ePrt  (58) 

s^Pi 


Since  the  Laplace  transform  of  ep'  r  is  — 2^,  y(t)  can  be 
transformed  into  the  Laplace  domain,  resulting  in  the  follow¬ 
ing  rational  form  of  the  transfer  function: 


7,  2  _ AC  : 

/C02  cCH3OH,ss^ell,ss 


GM3l(s)  —  ^02CcH3OH. 


a  — Pe  SS 


2  dM 


1  —  e 


oo  j 

Lat(.v)  =  Y/  lim  (x  -  P/)Y(s)) - - 

4-?s-+Pi  s-Pi 

1=1 


(59) 


■  [Pess  -  2  ePe^2F,(s)  +  2 P2(s)]  (50) 


nM  p  /)M 

r,  ,  ,  ^CH3OH  ,  ess  ;yCI[3()l[  „  ,  , 

Gm320s)  = - - V  K02*cell,ss  H - - F2(s) 

(51) 

with 

J(Pess/2)2  +  ((dM)2/D^ou)s 
Fi(s)  =  V  l  =  (52) 

sinh  f(Pe ss/2)2  +  ((u'm)2/£>cS3oh)5 


For  functions  with  an  infinite  number  of  singularities,  the 
higher  terms  represent  faster  dynamic  modes  (<  millisec¬ 
onds)  and  more  negative  singularities  P,.  These  become  ir¬ 
relevant  for  the  contemplated  order  of  magnitude  of  experi¬ 
mental  relaxation  times  (seconds-minutes).  Since  the  trans¬ 
fer  function  should  have  the  same  stationary  gain  as  Y(s ),  i.e. 
7(0)  according  to  the  final  value  theorem,  the  sum  in  Eq.  (59) 
is  convergent  for  s  =  0.  Therefore,  a  good  approximation  can 
be  achieved  by  truncating  the  series  after  the  time-relevant  m 
terms,  i.e.  time-relevant  m  singularities,  and  correcting  for 
the  stationary  gain: 


J(Pess/ 2)2  +  ((dM)2/D£C)s 
F2(s)  =  V  .  = 

tanh  J(Pe ss/2)2  +  ((dM)2  /  D^OH)s 


Y(s)  «  Yrat(s)  -  Yral( 0)  +  7(0)  (60) 

This  rationalising  procedure  is  applied  to  the  transcendental 
parts  Pi  (Eq.  (52))  and  P2  (Eq.  (53))  of  the  transfer  functions, 
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time,  t/[s] 


Fig.  15.  Comparison  of  voltage  responses  of  the  2x  +  m -model  (solid)  and 
the  2x-model  (dashed)  to  a  current  step  from  520  to  260  A  m~2  with  taio  = 
10-6  mol  m-2  s-1 ,  rA20  =  5  x  10-8  molm-2  s_1. 

Gmi  i(^')-Gm32(^'): 


time  t/[s] 

Fig.  16.  Concentrations  and  overpotential  responses  of  the  nonlinear  2x  + 
m- model  to  a  current  step  from  520  to  260 Am-2.  From  top  to  bottom: 
Methanol  concentration  and  surface  coverage  in  the  anode  catalyst  layer, 
anodic  overpotential,  methanol  flux  at  the  interface  membrane/cathodic  cat¬ 
alyst  layer  and  methanol  concentration  in  the  membrane  (10th  discretisation 
element  of  50  elements  in  total),  cathodic  overpotential. 


lim(— l)*+l 
Pi 


[Pe  ss/2)2  +  k2jt2\ 

(dM)2/D£fQH  J 


(61) 


F2(s)  F2,rat 


(Pess/2)2 +  A:27T2\  1 

(^M)2/^cS3 oh  J  ^ 


(62) 


p2(s)  does  not  fulfill  the  condition  in  Eq.  (55),  but  /’2,rat  is 
still  an  adequate  fit  to  the  original  function,  as  can  be  shown. 

The  new  transfer  functions  are  inserted  into  the  block  di¬ 
agram  (Fig.  9),  and  the  cell  voltage  response  to  the  exper¬ 
imental  current  step  is  simulated.  The  reaction  parameters 
are  identical  to  these  of  the  2x-model,  as  well  as  most  of  the 
other  parameters  (values  see  Nomenclature).  The  results  are 
displayed  in  Fig.  15  in  comparison  to  the  voltage  response 
of  the  2x-model.  As  can  be  seen,  the  relaxation  time  is  iden¬ 
tical  to  that  of  the  2x-model.  Observed  differences  are  a  2  s 
delay  time  until  the  maximum  voltage  is  reached  (discussion 
see  next  paragraph),  and  a  lower  steady  state  of  the  2x  +  m- 
model.  The  absolute  overshoot  is  slightly  larger  than  for  the 
2x-model,  but  it  is  still  small  compared  to  the  experimental 
voltage  response,  whose  overshoot  is  10  times  larger.  Hence, 
for  the  given  parameter  set  and  with  regard  to  the  experi¬ 
ment,  the  results  of  the  2x-model  and  the  2x  +  m -model  are 
qualitatively  equal. 

Observing  concentration  and  overpotential  changes  of  the 
numerically  evaluated  nonlinear  2x  +  m-model  in  Fig.  16 
(see  Section  3.2.6),  we  see  the  same  effect  as  in  the  2x- 
model:  methanol  concentration  dominates  the  dynamics  of 
the  cathodic  overpotential,  and  thus  the  cell  voltage.  Even 
after  50  s,  the  COx  coverage  is  not  in  equilibrium. 


The  dynamic  behaviour  of  the  flux  at  the  interface  mem¬ 
brane/cathode  shows  the  effect  of  the  introduced  osmotic  drag 
term:  methanol  flux  instantaneously  decreases  at  t  =  0  due 
to  the  reduced  current  density,  hence  due  to  lower  osmotic 
drag,  while  subsequently  the  anodic  methanol  concentration 
increases,  causing  a  rise  in  the  flux.  This  non-minimal  phase 
behaviour  leads  to  an  overshooting  of  the  cathodic  overpo¬ 
tential  response  and,  as  a  consequence,  the  above  presented 
overall  cell  voltage  response  (see  Fig.  15)  also  overshoots. 
Just  like  in  the  2x-model,  the  anodic  overpotential  is  quite  in¬ 
sensitive  to  changes  in  the  methanol  feed  concentration  and 
surface  coverage. 

3.2.6.  Discussion  on  validity  of  transfer  function  models 

The  last  two  sections  presented  dynamic  responses  of  the 
2x-model  and  the  2x  +  m-model.  The  responses  were  ob¬ 
tained  by  linearising  the  rigorous  model  equations  around  the 
steady  state  and  calculating  the  voltage  transient  to  a  current 
step,  using  Faplace  transformation.  As  can  be  seen  from  the 
results,  the  models  can  describe  the  observed  experimental 
behaviour  qualitatively,  but  not  quantitatively. 

Between  the  current  densities  of  520  and  260  A  m~2,  the 
experimental  polarisation  curve  is  approximately  linear  (see 
Fig.  2),  and  hence  the  system  could  be  regarded  as  a  linear 
one.  This  will  be  shown  in  the  following. 

For  comparison  with  the  linearised  transfer  function  mod¬ 
els,  the  nonlinear  system  answers  to  the  2x-model  and  the 
2x  +  m-model  have  been  calculated  by  numerical  solving  of 
the  respective  set  of  nonlinear  differential  equations  (see  also 
Figs.  14  and  16).  Both,  the  nonlinear  and  linear  system  an¬ 
swers  are  presented  in  Fig.  17  (top)  for  comparison.  As  can 
be  seen  for  both,  the  2x-model  and  the  2x  +  m- model,  volt¬ 
age  responses  of  the  linearised  and  nonlinear  systems  have 
the  same  shape  except  that  the  nonlinear  answer  is  shifted 
to  higher  voltages.  This  shift  is  due  to  the  deviation  of  the 
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Fig.  17.  Comparison  of  linear  models  with  non-linear  models.  (Top)  Voltage 
response  of  non-linear  models  ( 2x :  solid  black;  2x  +  m :  solid  grey)  and 
transfer  function  models  ( 2x :  dashed  black;  2x  +  m:  dashed  grey).  (Bottom) 
Relative  error  of  linear  compared  to  non-linear  responses  of  the  2x-model 
(black)  and  the  2x  +  m -model  (grey). 

polarisation  curve  from  the  tangent  of  the  polarisation  curve 
at  the  initial  current  density.  A  quantitative  analysis  of  the 
deviation  is  given  in  Fig.  17  (bottom),  where  the  relative  er¬ 
rors  of  the  linear  models  vs.  the  nonlinear  models  are  plotted. 
It  can  be  seen  for  the  2x-model  that  the  relative  error  in  the 
steady  state  is  quite  small  (2.4%).  Even  less,  approximately 
0.02%,  is  the  change  of  this  error  versus  time.  Therefore,  it 
can  be  concluded  that  the  linearised  2x-model  is  equivalent 
to  the  nonlinear  2x-model  for  the  given  current  step.  This 
may  change  with  increasing  current  step  size  where  the  devi¬ 
ation  from  the  polarisation  curve  will  increase.  The  relative 
error  of  the  linearised  2x  +  m -model  is  only  slightly  higher 
(^2.4%),  and  the  variance  of  the  error  is  in  the  same  order  as 
for  the  2x-model.  Hence,  also  for  the  2x  +  m-model,  it  can  be 
concluded  that  it  is  equivalent  to  its  non-linear  counterpart. 


Gc(s)  = 


Fig.  IB.  Voltage  response  of  the  2.r-model  (solid)  and  4.r-model  (dashed) 
to  a  current  step  from  520  to  260  A  m~2  with  rA \ o  =  10  6  mol  m  -  s-1, 
fA20  =  5  x  10-8  mol  m~2  s-1 . 

neously,  i.e.  in  a  quasi-steady  state.  This  assumption  seems 
justified,  since  in  the  experiments  the  Ohmic  part  of  the  volt¬ 
age  change  ( RM  ■  A/cen  =  3  mV)  is  only  a  small  part  of  the 
total  steady-state  voltage  change  (=52  mV)  for  the  standard 
current  density  step,  i.e.  from  520  to  260  Am'  2. 

Nonetheless,  if  the  given  set  of  equations  (Eqs.  (3)— ( 12)) 
and  parameters  (see  Nomenclature)  are  reliable,  a  model  in¬ 
cluding  the  dynamics  of  the  charge  balances  should  deviate 
only  marginally  from  the  2x-model.  In  the  following,  it  will 
be  shown  how  the  system  dynamics  changes  when  the  dy¬ 
namics  of  charge  balances  is  not  neglected.  Additionally  to 
the  two  state  variables  t'cH  oh  and  ^COx  of  the  2x-model, 
the  state  variables  i},\  and  >/c,  are  described  by  dynamic  bal¬ 
ances  (charge  balances,  Eqs.  (5)  and  (7)).  The  model  now 
consists  of  four  state  variables,  subsequently  it  is  denoted 
as  4x-model.  The  total  transfer  function  is  modified  by  ex¬ 
changing  the  transfer  function  elements  Gae  and  Gq  in  the 
blockdiagram  (see  Fig.  9): 

_ 1/Ccc _ \  T 

s—6/ Ccc F2roo(ac—1)/ R/T exp((ac— l)F/it/7)jc,ss)  | 

-6  F/Ccc  I  ^ 

s-6/CccF2rco(ac-l)/R/Texp((ac-l)F/R/Tr,c,ss)  ) 


GaeCv)  = 


/  -4FrA10exP(o!A^A,ss/«/r)/cf;ef3OHa-ecox,ss)/(i-ecox)/cAC  \ 

5-aA/J?/rexp(aAF),A,ss/t?/r)/CAC(-4F2rA1o^g3OHss/4c‘)3OH(l-ecOx,ss)/(l-0cox)-2f2''A2O^COx,ss/^ox) 

(4FrAloexp(aA^A,ss/t?/r)CAC30HiSs/cgf30H/(l-e«fx)-2FrA2oexp(«AF)7A,ss/fi/r)/e«fx)/CAC 

i-aA/R/rexp(«AF!7A,ss/t?/7’)/CAC(-4F2rA,0cAg3OHss/4c‘l3OH(l-ecOx,ss)/a-^Ox>-2f2''A20ecOx,ss/eg)x) 

_ 1/CAC _ 

\  i-aA/R/rexp(ffAF!7A,ss/R/7’)/CAC(-4F2rA|ocAC3OH  ss/c;:ef3OH(1_ecOXiSs)/(1_^f)[)_2f2,.A2oecOxiSs/0refx)  J 


(64) 


5.2.7.  Reaction-adsorption-charge  model  with  four  state 
variables  (4x-model) 

In  all  models  presented  in  the  previous  sections,  it  was  as¬ 
sumed  that  the  time  constants  of  the  charge  balances  are  much 
smaller  than  the  time  constants  of  the  mass  balances,  which 
means  that  the  charge  balances  are  in  equilibrium  instanta- 


Using  the  same  set  of  parameters  as  in  the  2x-model,  the  dy¬ 
namic  voltage  response  is  calculated  for  the  current  density 
step  from  520  to  260  Am-2.  The  voltage  response  of  both, 
2x-  and  4x-model,  is  shown  in  Fig.  18.  As  can  be  seen,  the 
time  for  reaching  the  new  steady  state  is  the  same  in  both 
models.  Five  seconds  after  the  current  step  both  answers  are 
identical.  But,  as  expected,  there  is  a  slight  discrepancy  in 
the  domain  of  the  small  time  constants,  i.e.  in  the  first  3  s.  In 
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the  4x-model,  the  voltage  cannot  step-up  immediately  to  the 
new  level  at  t  — »•  0+,  as  it  is  found  for  the  2x-model,  since  the 
double  layers  first  have  to  be  discharged.  This  discharge  takes 
place  within  three  seconds,  after  which  the  voltage  reaches  its 
maximum.  In  the  subsequent  time  only  small  voltage  changes 
due  to  concentration  changes  are  observed,  compared  to  the 
initial  large  voltage  changes  due  to  the  current  step,  and  the 
voltage  decreases  slowly.  As  a  result,  comparing  the  abso¬ 
lute  overshoots  of  the  2x-model  and  4x-model,  the  overshoot 
of  the  Ax- model  is  15%  smaller  than  that  of  the  2x-model. 
Nonetheless,  one  has  to  take  into  account  that  the  voltage 
response  will  always  show  this  characteristic  rounding  off 
of  the  overshoot  within  the  first  few  seconds  (time  depends 
on  the  double  layer  capacities),  but  the  4x-model  will  not 
contribute  other  dynamic  effects  to  the  system.  Finally,  the 
previously  postulated  assumption  that  the  4x-model  deviates 
only  marginally  from  the  2x-model  proved  correct,  the  set  of 
equations  showed  the  predicted  behaviour,  and  for  the  sake  of 
model  reduction,  the  models  with  quasi-steady  state  charge 
balances  may  be  used  for  modelling  further  on. 


4.  Conclusions 

This  work  presented  investigations  on  the  dynamic  be¬ 
haviour  of  the  DMFC.  Experimental  cell  voltage  responses 
to  load  step  changes  were  investigated  based  on  a  transfer 
function  approach.  The  here  presented  models  hold  for  anal¬ 
ysis  of  a  specific  operating  regime  (low  anodic  residence  time 
(~1  s),  low  current  densities),  in  which  dynamic  processes 
occurring  inside  the  MEA  are  minimally  influenced  by  mass 
transport  phenomena  (CO2  bubble  generation,  concentration 
distribution  in  anode  compartment). 

Applying  cell  current  steps  to  the  DMFC,  the  voltage  re¬ 
sponse  in  most  cases  showed  significant  overshooting  be¬ 
haviour.  The  height  of  overshooting  strongly  varied  depend¬ 
ing  on  the  position  of  the  investigated  steady  state  of  the 
DMFC  at  the  polarisation  curve,  hence  on  the  dominating 
physico-chemical  phenomena  in  different  operating  regimes. 
Switching  to  OCV  caused  overshoots  of  80%.  Understanding 
such  dynamic  effects  is  not  only  essential  for  DMFC  opera¬ 
tion,  but  also  for  understanding  the  processes  occurring  inside 
the  cell. 

Linear  system  analysis  was  applied  to  a  set  of  dynamic 
models  with  different  complexity.  A  criterion  for  cell  voltage 
overshooting  was  developed,  with  which  parameter  domains 
for  overshooting  were  analytically  determined.  Although 
even  a  model  with  only  one  state  variable  (i.e.  methanol 
concentration  in  the  anode  catalyst  layer)  is  able  to  generate 
cell  voltage  overshooting,  a  much  better  accordance  with  the 
experiments  was  achieved  when  introducing  a  consecutive 
reaction  mechanism  and  thus  working  with  two  state 
variables:  the  methanol  concentration  and  the  COx  surface 
coverage  in  the  anode  catalyst  layer.  The  implementation  of 
additional  dynamic  state  variables  (methanol  concentration 
in  the  membrane,  cathodic  and  anodic  overpotentials), 


for  the  given  parameter  set,  led  to  only  slightly  different 
results.  In  all  investigated  model  variants,  the  dominating 
physico-chemical  phenomena  were  the  same:  Transport 
of  methanol  through  the  membrane  caused  cathodic 
overpotential  overshooting  while  the  anodic  overpotential 
proved  to  be  rather  insensitive  to  methanol  concentration 
changes. 

Linearised  models  were  found  to  be  adequate  for  dynamic 
modelling  in  the  lower  pseudo-ohmic  regime  of  the  steady- 
state  polarisation  curve  at  high  anode  flow  rates,  hence  they 
also  could  be  suitable  for  controller  design  for  this  regime, 
which  will  be  the  scope  of  future  works.  Furthermore,  the 
causes  of  the  observed  overshoot  height,  which  could  not  yet 
be  reproduced  quantitatively  by  the  presented  models,  will 
be  investigated  (more  realistic  anode  and  cathode  kinetics). 

The  transfer  function  library  enabled  the  modular  combi¬ 
nation  of  descriptions  of  physico-chemical  phenomena  inside 
the  MEA  of  the  DMFC.  Further  investigations  will  be  con¬ 
ducted  by  enlarging  the  library  with  modules  describing  outer 
mass  transport.  Finally,  the  transfer  function  library  will  be 
used  for  EIS  simulations. 
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